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ABSTRACT 


The statistics of the fluctuations of the in-situ speed of sound 


in the upper ocean have been Studied by analyzing the instantaneous 





phase difference of the Output of two hydrophones Separated by one 
meter for sounds of frequency 15 to 151 kHz. The experiment was con- 
ducted at 11 ft in water of depth 60 ft in low sea states at night. 
Comparison of the speed calculated from the time averaged phase dif- 
ference, with the speed given by velocimeter or empirical se oa 
yielded differential speeds which deviate by 1 m/sec to 8 m/sec from 


the accepted values, for frequencies less than 100 kHz. Correlation 
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and spectral analysis of the sound phase and ocean height fluctu- 
ations has shown the Close relation between these two parameters. 
There is strong evidence of the presence and importance of bubbles 
in all of the results, particularly of a large population resonant 
in the frequency range 56.3 to 71.1 kHz (radius 50 to 60 microns). 
Evidence is presented to Suggest that bubbles appear at the surface 
during internal wave activity at lower depths and that for sound 
peumnlenio ee near the bubble resonances the sound phase is strongly 


helices by, and in phase with, the ocean wave height. 
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I. INTRODUCTION 


The objective of this research was to measure acoustic phase 
shifts over a one meter path in order to study the statistics of 
the variation of sound speed with frequency and time in the upper 
layer of the ocean. 

The sources of the dispersion of sound speed are bubbles of 
various radii. Motion of these bubbles and of other inhomogenieties 
cause the temporal fluctuations of phase. 

Study of the dispersive behavior is important because in the 
top tens of meters of the ocean, the speed of sound at military 
frequencies differs both from the value given by velocimeters 
(which operate in the megahertz range) and that given by the vari- 
ous accepted empirical formulas. Study of the fluctuations is 
necessary in order to reveal the oceanographic origins and de- 
pendencies of the temporal variations in sonar surface duct 
propagation. 

The bubble problem has been treated theoretically by Meyer 
and Skudrzyk (1950), Devin (1959) and Albers (1960). The experi- 
mental study of wind generated bubbles has been conducted in the 
laboratory by Glotov (1962) and at sea by Buxcey et al (1965) and 
Medwin (1970). Most recently Rautmann (1971) measured the magni- 
tude of dispersion and fluctuation of sound speed in the surface 
layer of the ocean. 

Change in the speed of sound with frequency in sea water 
is attributible in general either to molecular relax- 


ation or bubbles since no other frequency dependent mechanisms 
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are known to exist at the frequencies of this study. The magnitude 
of the molecular relaxation can be shown to be negligible as com- 
pared to the bubble effects (a frequency shift from 20 to 70 kHz 
results ina Ac of only approximately 0.3 m/sec). 

The research was conducted using the U.S. Naval Undersea 
Research and Development Center (NUC) Oceanographic Research Tower 


in Mission Bay, San Diego, as a test platforn. 








II. INSTRUMENTATION 


An aluminum pipe frame of height six feet and length twelve 
feet was used to mount the transducers and hydrophones (Figure 1). 
To prevent displacement of the hydrophones, the frame wires were 
put under tension by spring loaded mounts. To minimize reflections 
from the frame caused by side-lobes of the transducers, the frame 
was pulse-echo checked; various segments of the frame were covered 
with acoustic absorbent rubber (SOAB) so that the echo was down 
35 dB below the direct signal. 
| Opposite to the open end of the frame, at the center, and 
pointing in the X-direction, an USNRD F33 unidirectional trans- 
ducer was mounted. On the acoustic axis of this transducer, two 
Atlantic Research LC-10 hydrophones were placed. The first LC-10 
was mounted 82.3 cm from the source and the second was mounted 
81.3 cm beyond the first. 

The transmitting instrumentation (Figure 2) consisted of a 
General Radio Coherent Decade Frequency Synthesizer, type 1162-A, 
generating a stable sinusoidal voltage output of two volts rms 
and variable in frequency. The signal was amplified by a Hewlett 
Packard Power Amplifier, 467A, to 18 volts rms and then impressed 
across the transducer. 

The receiving instrumentation (Figure 2) consisted of LC-10 
#1 with a flat (-108 dB + 2 dB re 1 volt/microbar) receiving re- 
sponse from 15 kHz to 150 kHz, and LC-10 #2 with a flat (-110 dB 
+ 3 dB re 1 volt/microbar) receiving response from 15 kHz to 


150 kHz. The signals were then amplified 30 dB by an NUS 
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pre-amplifier model 2010-030 before passing through 150 feet of 
waterproof shielded cable. Both Signals were then bandpass filtered 
by a digitally tuned variable Krohn-Hite model 3322 Filter set with 
the center frequency equal to the transmitted frequency + 300 Hz. 
The filtered signals were then fed into a Dranetz model 305 Phase 
Meter. The phasemeter measured the phase angle between the two 

ac voltages at the same frequency and provided a dc voltage output 
proportional to the difference of the phase angles at 10 mV per 
degree. The dc voltage was Simultaneously FM recorded ona 
Precision Instrument model 6200, eight track magnetic tape re- 
corder at 3.75 1.p.S. and on a Brush model 220 strip chart 
recorder, at 10 mV per division and 1 mm per second. 

The ocean wave height was measured utilizing the installed NUC 
Baylor wave height gauge. The wave height gauge was mounted a 
distance of 6 feet from the acoustic source and at an angle of 
60 degrees from the acoustic axis. The dc voltage output was FM 
recorded with each phase record and simultaneously recorded on the 


strip chart recorder. 
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III. OCEANOGRAPHIC COMPONENTS 


Field measurements were similtaneously taken to determine the 
- fluctuations of temperature, salinity, water particle velocity and 
wave height to obtain a better understanding of the acoustic-ocean 
interaction (theses by CDR J. Gossner, LT M. Whittemore, 

LT W. Frigge and LT R. Krapole, March 1973). 

Unfortunately much of this information is unavailable at this 
time. However, the time-varying ocean wave height, (NUC Baylor 
Gauge), the sound velocimeter output, (Ramsey Probe), and tempera- 
ture (Ramsey Probe) were simultaneously recorded during the phase 
study. The average value of salinity (33.8 °/00) measured that 
afternoon was used for the C calulations on the assumption that 


changes in Salinity would have a minimal effect on C. 
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IV. EXPERIMENTAL SETUP 


The frequency of the continuous wave signal applied to the 
transducer was adjusted by means of the frequency synthesizer 
until a phase angle of approximately 20 degrees was obtained. 

This value was selected to prevent the possibility of the phase 
going negative and changing the integral number of wave lengths 
between the hydrophones. Since the phase meter compares the signal 
at input B with respect to the signal at input A, and cannot de- 
tect integral numbers of wave lengths difference between the two 
signals, the displayed phase determines the fractional wave length 
by which Signal B leads signal A. This means an integral number 
of wave lengths plus a fractional part fits into the separation 
distance between the two hydrophones. Increasing the frequency 
will change the phase by 360 degrees or step the integral number 
of wave lengths by one. 

Because it was desired to study the statistics of the speed of 
sound as a function of frequency over a range most likely to con- 
tain resonant bubbles, the frequency range from 15-150 kHz was 
examined. The 81.3 cm separation between the two hydrophones 
caused the number of wave lengths within the spacing to be incre- 
mented by one approximately every two kHz. 

In order to exclude possible errors in the results due to non- 
stationarity in the medium, it was necessary to examine the total 
frequency range in as short a period of time as possible. The 


The frequency range was examined in steps of four kHz, and the 
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time variable phase difference and wave height were recorded be- 
tween the hydrophones for four minutes at each frequency. 

All the readings were taken within three consecutive hours, 
from 1900 to 2200 hours on June 5, 1972 at an initial transducer 
depth of 3.4 meters, and final depth of 2.8 meters. The tide 
during this period was falling at a rate of 0.5 ft/hr, with low 


water at approximately 2330. 
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V. DATA REDUCTION 


Thirty-four four minute individual physical records make up 
the analog recorded sound phase and wave height data. One ap- 
proach used was to digitize the analog data and use a computer to 
carry out the statistical and spectral analysis. 

The recorded data were digitized using the NPS CDC 5000 and 
SDS 9300 computers. The initial sampling frequency was 50 Hz 
giving a Nyquist frequency of 25 Hz. The sampling interval was 
0.02 seconds. The 7 track digital tape was then processed to 
remove spurious signals (Figure 3), scaled, and conver ted using 
the IBM 360/67 from 7 track octal to 9 track hexidecimal for use 
by the IBM 360/67, (computer program 1). 

The samples generated, however, by uSing the 50 Hz sampling 
rate provided too much data to analyze to the degree of resolution 
required in any reasonable processing time on the IBM 360. It was 
then determined to reduce the rate to 6.25 Hz, a sampling interval 
of 0.16 seconds, and frequency resolution of 0.0244 Hz. This 
yielded 10 degreesof freedom. This also reduced the amount of each 
physical record that could be analyzed from 3 minutes 50 seconds 
to 3 minutes 40 seconds. 

A second approach was essential in order to obtain statistics 
of the average and standard deviation which would be a measure of 
the microstructure rather than of the microstructure plus the 
internal wave. To minimize the effect of internal waves each 


record was analyzed in 30.72 second increments. Each increment 
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consisted of 1536 samples for which instantaneous values of the 

velocity of sound were calculated and the means and standard de- 
viation were determined. Similarly the mean value and standard 

deviation of the phase were computed, (computer program 2). 

The speed of the coherent component was also calculated using 
the mean value of the phase and compared to the mean value of the 
speed. The point being made is that, if the mean of the entire 
physical record had been used, the large variance would have been 
relative to a mean value strongly influenced by the passage of 
any internal waves and not due to the microstructure of the medium 
as desired. The mean speed and the mean coherent speed of sound 
were then compared with the 30 second average values given by a 
Ramsey Velocimeter and those calculated by empirical formulas 
uSing the salinity, depth, and 30 second averages of the 
temperature. 

For selected frequency runs, the temporal autocorrelations and 
power spectra, in dB re degrees, or dB re cm, of the phase and wave 
height were Saameeed and plotted using a maximum time lag of LO .96 
seconds. The temporal cross-correlation was also computed. From 
this the cross spectrum, coherence, and cross spectral phase angle 
were computed, (computer programs 3 and 4). The selected fre- 
quencies for which these were calculated were determined by first 
calculating and plotting the dispersion curve (Figure 4), and then 


selecting those frequencies of interest from the results. 
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VI. COMPUTATIONAL METHOD 


The computational method for determination of the speed of 
sound using the phase difference was the same as that used by 
Rautmann (Dec. 1971) and is included here for information. 

Since a continuous wave method was used, the signal received 


by the first hydrophone can be written as: 
i sin (Wt) 


and the signal received by the second hydrophone is: 


b sin (Wt-kx) 


< 
i 


B 
where 
k = the wave number = ie oe 
Cc Cc 
f = the frequency of the applied signal 
C = the phase wave velocity 
x = the separation between the hydrophones. 


Since we compare the two signals with respect to phase by 


using the digital phase meter we have to find values for which 
kx = U 


where YW is the total phase difference between the two signals due 
to an integral number of wave lengths and a fraction of the wave 


length. This is written as: 


kx = (27)n + Ph n = an integer 
a3) 
kx = = = ; 
x Fox 201) nee go. 9 in radians 
Os 
C= 


” yn +B, 
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The final result is: 





where is in degrees. 
n+ @ 2 g 


The integer n can be pre-estimated by calculating C from Wilson's 
equation uSing the known values of the temperature and salinity. 
Therefore the instantaneous speed of sound may be obtained by 
uSing the measured frequency separation and whole number of wave- 
lengths at which a residual phase difference is noted. 
The value for the speed of sound of the coherent component 


is defined to be: 


Sees 
a =o. where (@) is the mean value of the phase. 
360 
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VII. DISPERSION THEORY 


The following equation relates bubble size to the speed of 
of sound and therefore shows the dispersion for a bubble popu- 


lation of unique radius R. 


zZ =a | 
1 By o U(R)dR QQ 


os 20 R- G) - 79) 5 
ered #] + Bo =e 
) ) 


nondispersive sound speed 


Cc 
oO 


G) = sound frequency 


Cu(R)dRJ= fractional air to water volume for the increment between 
R and R + dR 


& 


= resonant frequency of bubble of radius R 


damping constant of bubble of radius R. 


A typical dispersion (differential sound speed) curve is 
shown below. A graph of the variance of C due to a fluctuation 
of the number of bubbles for a single size bubble as studied by 
Professor P.C.C. Wang (personal communication) would be M-shaped. 

\ 
In addition, a constant number of bubbles of changing radius (say, 
due to changing pressure) would create a variance of C which 
should peak at the mean resonant frequency fo: as shown by the 
dashed curve. 

A mixed bubble population would have the same equation summed 
over all of the bubble radii. The various concentrations of 
bubbles of different resonant frequencies in the ocean would cause 
combinations of graphs such as shown on the following page. The 


Ws, 





Ceo 


VAR C 


where f = 
Oo 

vy = 

Vo = 

6 = 

Var C= 








ae em eee eee eee 


Fig. 5 


mean resonance frequency of bubble 


small bubble density 


large bubble density 


damping constant of bubble 


variance of speed of sound due to fluctuation of 
the number of bubbles 


Five distinctive regions of the simple case shown are identi- 


fied for future use: 


Region 
Region 


Region 
Region 


Region 


A: 
B: 


Ce 
D: 


low variance, small negative differential speed 
maximum variance, maximum negative differential 
speed 

large variance, minimum differential speed 
large variance, maximum positive differential 
speed | 


E: low variance, small positive differential speed. 
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differential sound speed and the variance of the sound speed are 
studied as functions of frequency in order to deduce the bubble 


concentrations of various sizes. 
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VIII. EXPERIMENTAL RESULTS 


The frequency run was conducted at a depth of 3.4 meters at 
1900 decreasing to 2.8 meters by 2200. The run covered a total 
frequency range of 15 kHz to 151 kHz in steps of approximately 
4 kHz. 

Continuous direct recordings of temperature and sound velocity 
were taken throughout the frequency run by a Ramsey Velocimeter. 
Because of the internal wave presence, Ramsey velocities and tem- 
peratures were averaged over 31 second segments for comparison 
with the speed of sound measured by the phase method in the same 
interval. During that same interval, calculations of the empirical 
values of the speed of sound were done uSing both the generally 
accepted Wilson's equation {1] and the more recently proposed 
Leroy's equation [2]. Both methods yielded results within 0.1 
m/sec over the range of physical parameters of the experiment 
(see Table I). 

Wilson's equation was used in the following form: 

c = 1449.14000 + 4.57210(T) - 0.04453T7 ~ 0.00026T° 


+ 0.000017! 


+ 1.39799(S-35.0) + 0.00169(S-35.0)* 
+ 0.06240 + (S-35.0)(-0.01244T). 
Leroy's equation was used in the following forn: 
C = 1492.9000 + 3(T-10.0) - 0.00600 (T-10.0)* - 0.04000 
(T-18.0)7 + 1,20000(S-35.0) -0.01000(T-18.0)(S-35.0) 


2 


+ Z/61.0 + 0.1RHO” + 0 .00020RHO* (T-18 .0) * + 0.10000RHO 


(8/90.0) + 2.0/(10.0) /(T)(T-10.0)4 


Ze 





where in 


RHO 


temperature in 


the equations: 


salinity in parts per thousand 


depth in meters 


latitude 


z/1000. 


degrees centigrade 


speed of sound in meters per second 


The Leroy equation bubble free values were used as the refer- 


ence for the differential speed of sound curve (Figure 4). The 


multi-peaked differential sound speed curve suggests several dif- 


ferent bubble concentrations of different resonant frequencies. 


Therefore, the curve was studied and the following frequencies 


were selected for analysis: 


Frequenc 


14.79 
18.78 
26',17 
41.32 
56525 


63.79 


71.13 


89.90 


112.52 
123.72 
136.80 
146.10 


TABLE II: Selected Frequencies 


kHz Reason for Selection 


minimum speed, 


maximum speed, 


possible zero 
maximum speed 


possibly near 
frequency 


maximum speed 
minimum speed 
minimum speed 
minimum speed 
minimum speed 
minimum speed 


minimum speed 
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lowest frequency 
low frequency 


dispersion frequency 


the zero dispersion 


at higher frequency 
at higher frequency 
at higher frequency 
at higher frequency 


at higher frequency 
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A. DIFFERENTIAL SPEED OF SOUND AND VARIANCE OF SOUND SPEED AS A 
FUNCTION OF FREQUENCY 


The variation of the speed of sound with frequency was determined 
by subtracting the values calculated by Leroy's equation (a function 
of temperature, salinity, depth and latitude) from the experimental 
values. 

At each frequency a pre-run electrical calibration was conducted 
to determine the phase error correction due to differential phase 
shifts in the amplifier and band-pass filters. This correction was 
then algebraically subtracted from the phase measured during the ex- 
periment to obtain the corrected phase difference between the sig- 
nals at the hydrophones. 

There were two other sources of possible corrections that were 
investigated. The first was motion of the medium, such as an onshore 
current. In a previous experiment, it was determined that there was 
no de flow at this point. However, any correction would have been 
small and would affect only the absolute values and not the relative 
values. This correction would be independent of frequency. 

The second was measurement error in determining the exact dis- 
tance between the acoustic centers of the two hydrophones. A mea- 
surement error of 2.2 mm would yield a 4 m/sec change in the speed 
of sound at a frequency of 136.80 kHz. This measurement error 
would affect only the absolute values of the speed of sound. How- 
ever, the relative values of the differential speed of sound are 
of interest, not the absolute values. Therefore this correction 
is not relevant. The dashed line on Figure 4 illustrates the 
effect on a 2.2 mm error in the measurement of the distance be- 


tween the hydrophones. It is believed that this probably is close 


to the correct absolute value as the speed of sound correction 
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would be small at high frequency. In theory, the bubble effect 
becomes negligible as the frequency increases (section VII). As 

the frequency increases, the speed of sound calculated by the phase 
method asymptotically approaches the values predicted by the Wilson 
and Leroy equations. This is also the reason the Ramsey velocimeter 
values agree with the predicted values. The Ramsey Velocimeter 
operates in the megahertz range, well above any bubble effects. 

The speed of sound values were calculated, and the differential 
sound speed determined, by subtracting the Leroy's equation values 
from the phase calculated values. 

The values obtained from Leroy's equation were in agreement 
within + 0.1 m/sec of those obtained from Wilson's equation. The 
Ramsey velocimeter values were within + 1 m/sec of both the Leroy 
and Wilson values. The calculated values from the four methods 
are included as Table I. 

The values are plotted in Figure 4. What are of interest are 
identifiable bubble concentrations of various size bubbles. The 
graph is similar to the one presented by Rautmann, based on mea- 
surements at the same location in October 21, 1971. Troughs and 
peaks of speed, presumably bounding a relative zero differential 
speed and loosely corresponding to the regional identification 
BCD defined in Figure 5, appear to be present at frequency bands 
roughly 15-19 kHz, 22-40 kHz, 50-65 kHz, 70-80 kHz, 100-120 kH2. 
This tenative description will be confronted with other statis- 
tical evidence, in sections VIII A and B. 

The twelve selected frequencies of Table II. will now be 


considered in detail: 


26 





The frequency 14790 Hz shows a corrected differential sound 
speed of -12 m/sec. This suggests that this frequency is below 
the resonant frequency of a large concentration of bubbles. Theory 
predicts a maximum standard deviation in this frequency range 
(Region A) caused by bubbles resonant at any higher frequency. The 
large standard deviation calculated was 1.57 m/sec therefore again 
suggesting that this frequency is in the range f < fo: 

The frequency 18777 Hz shows a corrected differential speed of 
-2.4 m/sec. If the (-) sign is correct, this would suggest 
region A of the bulk of bubble population but in region C of a 
large concentration of resonant bubbles at an f less than 18.9 kHz. 
If this is true, a batch of large bubbles (radius greater than 
200 microns has been identified, which could only have come from 
the bottom, or which perhaps represent biological entities. The 
large standard deviation of 0.92 m/sec may show that this fluctu- 
ation in fy is indeed taking place. 

At 26170 Hz, the corrected differential sound speed is -5.7 
m/sec. This may imply region B of a large population of bubbles. 
The standard deviation is large at 0.84 m/sec. 

At 41320 Hz the corrected differential sound speed is +2.6 
m/sec. The standard deviation is 0.58 m/sec. This value of the 
standard deviation implies this frequency is above the resonant 
frequency of a major bubble concentration. 

The 56245 Hz record shows strong evidence of a Slowly changing 
speed of sound, such as for an internal wave. The internal wave 
inferred from phase information reached a phase trough approxi- 


mately 40 seconds into the 4 minute record. This resulted in 
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measured instantaneous phase angles on the order of + 13 degrees, 
at the beginning of the record, to + 25 degrees at the passage of 
the peak of the wave. During this period, the values of the 
instantaneous sound speed measured by our phase difference techni- 
que decreased 1.7 m/sec while the Ramsey Velocimeter showed no 
appreciable change. That this is not clearly reflected in Table I 
is due to the 31 second intervals, chosen in the attempt to eli- 
Minate these effects, over which the means and statistics were 
computed. More will be said about this apparent discrepancy later. 
The results here eee indicate evidence of a resonant bubble 
population. The corrected differential sound speed is + 2 m/sec 
and the standard deviation is large at 0.55 m/sec. This large 
value for the standard deviation may mean this frequency is in 

or near an fo: As in the 18777 Hz record a small fluctuation in 
fy would cause this large OO: 

The 63790 Hz record has the maximum corrected differential 
sound speed value of + 3.8 m/sec. The standard deviation is also 
large, 0.55 m/sec. This may mean that a frequency close to the 
resonant frequency of a large bubble population. 

At 71130 Hz the record shows the largest local value in the 
standard deviation of the speed, 0.71 m/sec, (Figure 6). If the 
variance is due to change of number of bubbles it would represent 
a peak in the standard deviation at £ N1+6 or f N1-6. However , 
if near a concentration of resonant bubbles, whose frequency is 
changing in Region C, the fluctuations of the surface wave and 
upper ocean particle velocity may also cause a large variance. It 


remains then to investigate the correlation functions to interpret 


28 





(a. | 


STD. DEV. 
M/SEC 


O 
O 


O 
iN 


S 
ND 





20 60 100 140 


mee OUaNE YG) 
2 ANE TIME 1900 


Figme 


29 





this record. In any event, since 6 ™ 0.1, the peak standard 
deviation at a frequency is within 5% of the resonance of a large 
bubble population. 

At 89900 Hz, the corrected differential speed of sound is down 
to -1.4 m/sec and a small standard deviation of 0.28 m/sec. This 
suggests a resonance above 89.9 kHz. The remaining records at 
112519 Hz, 123720 Hz, 136800 Hz, and 146100 Hz all have small 
standard deviations and small corrected differential sound speed 
values. This would indicate that these frequencies are above the 
region of any major bubble populations. As the frequency in- 
creases, in theory, the sound speed should approach the empirically 
predicted or velocimeter values. That the speeds are now rea- 
sonably constant in the frequency region 130-150 Pemetedes ce that 
the sound speed is approaching this asymptote, perhaps with a small 
error in the 4 m/sec correction. This suggestion, in the case of 
the sound speed, is supported by the now slowly increasing value 
of the variance of the phase as the frequency increases from ap- 
proximately 100 kHz to 150 kHz. An increaSe proportional to ¢° 
is predicted by Chernov (equation 141, page 75, Ref. 5), for a 
medium without bubbles. At frequencies above 100 kHz, it appears 
that movement of the sea surface has little effect on the sound 
phase except insofar as the wave induced well correlated particle 
velocities underwater move temperature inhomogeneities through 


the sound path. 
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B. OTHER STATISTICAL ANALYSES 
1. Probability Density Functions 

Computer program 2 was used to compute the instantaneous 
values of the speed of sound and to form histograms for each of 
the 30.72 second intervals. This allowed interpretation of the 
speed of sound distribution functions. The predominant pattern 
of the envelopes followed a Gaussian distribution for the lower 
of the frequencies. This was expected because in the measurements 
many weak fluctuations in temperature and salinity would cause 
small random fluctuation. 

In studying the distributions of the selected frequencies, 
a change is observed at frequencies where the probability of re- 
sonant bubble populations were noted. At 14790 Hz the distribution 
is well defined and sharply Gaussian. At 18777 Hz, 26170 Hz, and 
41130 Hz, the shape is still that of a Gaussian distribution but 
now flatter with wider skirts. 

At 56145 Hz, the frequency where clear effects of bubbles 
were noted, the PDF begins to change the shape of its envelope. 
The envelope now shows two peaks in the distribution of velocities. 
Because of the presence of the internal wave it is difficult to 
determine whether this is predominantly due to the wave effects or 
the proximity of the resonant bubble population. 

The envelope at 63790 Hz also shows this double peaked 
behavior. This may again be due to the sound frequency being in 
the proximity of the resonance frequency of a major population of 
bubbles. The passage of these bubbles back and forth through the 


acoustic path and their size modulation due to the changes in the 
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hydrostatic pressure as the wave passes over are suggested as 
the cause of the non-Gaussian PDF. The PDF's at 71130 Hz also 
display this same behavior. 

The PDF behavior above these frequencies 1s not clear. At 
89900 Hz a double peaked envelope appears in 4 of the 7 segments 
analyzed. At 112519 Hz the distribution again appears GausSian. 
In the remaining 3 records studied, the envelope was not con- 
sistant, one segment of a record being Gaussian and the next 
double peaked. It is less likely that bubbles are important here, 
which suggests that the velocity distributions at higher sound 
frequencies are the result of influence by internal waves. Because 
of the relatively diffuse nature of the data from the 84 PDF graphs, 
they have not been reproduced here. 

2. Correlation Analysis 

The temporal autocorrelation function is the normalized 

temporal autocovariance function, normalized by dividing the 


autocovariance function by the variance, its maximum value at 


T = 0. The temporal autocorrelation function is then: 
(t)x(t+T) 
R(T) = (x t = t ) 
0 
x 


The plots of the autocorrelation functions of the phase and wave 
height are included as Figures 7 - 18. 

Similiarly the temporal crosscorrelation function is the 
normalized temporal crosscovariance function, normalized by 
dividing the crosscovariance function by the square root of the 
product of the maximum values of the two signals. The temporal 
crosscorrelation function is then: 


R(T) = Cx(t)y (t+ . {ont )e 
low oF! : 
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The plots of the crosscorrelation function of the phase with the 
wave height are i ana as Figures 19 - 30. 

At 14790 Hz (Figure 7) the phase and wave height remain 
correlated only for approximately the first 3 seconds. The cross- 
correlation (Figure 19) of the phase with the wave height has a 
peak to trough value of 0.45, and is good over only a time of 
approximately 25 seconds. The periodicity of the crosscorrelation, 
approximately 16 seconds, corresponds to a surface swell and 
therefore a particle motion through the sound field. The cross- 
correlation does not peak at fT = 0, which suggests that there is 
a time lag between those two effects. This time lag could identify 
transport of bubbles into and out of the sound path and represent 
the time difference between peak wave height at the wave height 
probe 6 ft away and the maximum particle velocity of the sound 
path. 

At 18777 Hz (Figure 8) the phase and wave height remain 
relatively in step. The crosscorrelation (Figure 20) is good for 
approximately 16 seconds and has a peak to trough value of 0.4. 
This again suggests a resonant bubble population as implied in the 
previous discussion of this frequency in Section VIII A. 

The 26170 Hz (Figure 9) correlation functions show a phase 
periodicity on the order of 12 seconds, while the surface wave 
height periodicity was again 16 seconds. The crosscorrelation 
(Figure 21) function is not a maximum at T = O but has a high peak 
to trough value of approximately 0.6. 

The 41320 Hz (Figure 10) correlation functions of the 


phase and wave height are uncorrelated indicating that this 
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frequency was not at a resonant frequency of a bubble population. 
The crosscorrelation (Figure 22) has a peak to trough value of 
only 0.16. There is evidence of good crosscorrelation with an 
internal wave of periodicity perhaps 200 seconds. 

At 56245 Hz (Figure 11) the correlation functions are in 
step for a time lag up to approximately 28 seconds. Note the peak 
of crosscorrelation at zero time lag. The crosscorrelation 
(Figure 23) function has a peak to trough value of 0.6 at the 
swell average periodicity and also shows the effect of an internal 
wave. This is strong evidence of a resonant bubble population 
superimposed on internal wave effects. It was observed previously 
that the Ramsey Velocimeter showed no evidence of an ee wave 
at this time. In fact, the phase measurement 1S so close to the 
Bom falce), and in the mixed layer, that large temperature changes 
are not expected. The contradiction can be resolved if it is 
speculated that bubbles of resonant frequency 56245 Hz were gene- 
rated in step with the internal wave activity in the thermocline. 

The 63790 Hz (Figure 12) record shows a region of bubble 
concentration is near. The peaks of the phase correlation are 
mostly in step with the wave height; only the magnitude here 
differs. However, the phase correlation is not as narrow banded 
as the wave height and shows some evidence of an internal wave. 
The wave height displays the same narrow band and magnitude as the 
previous records. The peak to trough crosscorrelation (Figure 24) 
value occurs at zero time lag and has now increased to 0.8, the 
second largest value of the experiment. The double peak in the 


phase record near 30 seconds shows up here as in the 41320 Hz 
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functions, and is apparently a component of the sea surface 
height spectrum. This is verified by observation of the spectrum 
(Figure 37). 

At 71130 Hz (Figure 13) it is found that the temporal 
correlations of the phase and wave height are of equal narrow 
band and have the same values at the same lag times. The cross- 
correlation function (Figure 25) has a peak to trough value of 1 
and is in step at T = 0. This is the strongest evidence that the 
sound phase is changing in step with the ocean hydrostatic pres- 
sure. If the analysis is correct, this means that the phase is 
increasing and decreaSing as the wave height increases and de- 
creases. That is to say, this frequency is at the resonant 
frequency of a major population of bubbles and the speed is moving 
up and down the steep portion of the slope of the dispersion curve 
(Region C) for these prominent bubbles at this frequency. 

The correlation functions at 89900 Hz (Figure 14) show 
that the phase at this frequency is poorly correlated with the 
wave height. The crosscorrelation function (Figure 26) shows 
almost no correlation with wave height. There is some evidence 
of long period correlation with an internal wave. This means that 
the sound phase is almost undisturbed by the surface wave or the 
orbital particle motion. This implies this frequency is not near 
any major resonant bubble population frequencies. Note also the 
sound phase and surface wave height are no longer in step T = O. 

The wave height temporal Se eee function at 112519ehz 
(Figure 15) is particularly simple and that of a narrow band swell. 
This 1S a good example of a simple damped cosine auto-correlation 


function. The phase correlation is poor but shows small peaks 
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at 7, 12, 19, 28, 32, and possibly 40 seconds. The cross- 
correlation function (Figure 27) shows very weak correlation 
between sound phase and wave height. The peak correlation does 
not occur at T = 0. 

At 123720 Hz (Figure 16) the phase remains correlated with 
itself longer than the wave height and shows some contribution due 
to the wave height periodicity at 16 seconds. Additional phase 
correlation peaks show up at 26 and again possibly 40 seconds. 

The crosscorrelation (Figure 28) has a relatively strong cor- 
relation for approximately 25 seconds. The peak to trough value 
is 0.65 but it now shows a full 180 degree phase shift between 
sound phase and wave height. This means that instead of phase 
increasing with wave height, it decreases. This is interpreted 
to mean Region E for all bubbles, that is, there are no signifi- 
cant bubbles resonant above this frequency. 

In the remaining two records, 136800 Hz (Figure 17) and 
146100 Hz, (Figure 18) which were studied for the higher fre- 
quency effects, the phase again remains correlated with itself 
longer than the wave height, that is, it is relatively uninfluenced 
by wave height. The crosscorrelation has a weak peak to trough 
value of approximately 0.3. This means again that the phase is 
relatively undisturbed by the surface waves or that the sound 
speed is undisturbed by the particle motion. The 136800 Hz phase 
PeanSiation shows small peaks at 12 and 18 seconds. The 146100 Hz 
phase correlation shows small peaks at 10 and again at 18 seconds 
of lag. 

The peak-to-trough range of crosscorrelation values 


(Figures 29, 30) at the swell frequency, constitutes an important 
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parameter showing the ability of the sound phase to follow the 
surface wave height (or orbital particle velocity). A curve of 
this variation with frequency is shown in Figure 3l. 
3. Spectral Analysis 

The power auto-spectral density of the phase and wave 
height records gives the distribution, in the ene NsTICS domain, of 
average power of the fluctuations vs. frequency. The method used 
was to compute the autocovariance function and take an FFT and 
applying a Hamming window. This was done uSing computer program 
(3) and yielded plots of phase angle squared and wave height squared 
vs. frequency in Hz. The frequency resolution is 0.024 Hz. To 
improve the pictoral representation in amplitude, 10 109, 9 of the 
ordinate values was taken so that the power spectrum 1s a power 
spectrum ewe im dB re 1 a 

For all records of the wave height, there is a distinct 
peak at 0.06 Hz due to the observed ocean swell component, and 
some smaller peaks up to 0.4 Hz. The peaks up to approximately 
0.4 Hz could be due to a smaller wind generated wave superimposed 
on the larger swell component. The constancy of the spectral 
characteristics below 0.4 Hz, implies that over the 3 hour duration 
of the experiment, the ocean was not changing significantly. How- 
ever, the higher frequencies (above 0.4 Hz) with relatively con- 


stant levels are apparently some form of generated noise. 


An error appears to have slipped into the ordinate of Figures 
36, 40, 41, 42. The wave height power levels can be corrected 
by subtracting 26 dB from these ordinates. The visually 


observed rms wave height of approximately 2 ft, agrees with 
this correction. 
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The phase spectral records are not as consistent as the 
wave height records. The surface wave component contribution at 
0.06 Hz is readily identifiable in all except the 41.3 kHz record. 
Beyond that, any consistency is not so clear. The phaSe plots are 
in relative dB only. 

The peak at 0.4 that was observed in all the wave height 
spectrums was only observed in the lower three frequency phase 
records analyzed. The 0.4 Hz peak was again probably a surface 
wave component. In the vicinity of bubble population effects, 
this frequency could not be identified. 

The coherence and cross spectral phase angle is computed 
(computer program 4) by taking the Fourier Transform of the cross- 
covariance function and finding the cross-spectral density using 
a Parzen lag window. This yields a real part; the co-spectrun, 
and an imaginary part; the quad-spectrum. The magnitude is then 
the square root of the sum of the squares of the co and quad- 
Spectra. The square of the magnitude is then divided by the 
product of the individual spectral densities to yield the coherence 
function. The values so obtained range from -l to +l. The cross 
Spectral phase angle is found by taking the arctangent of the 
quotient of the quad-spectrum with the co-spectrum. The maximum 
values are then + 180 degrees. Interpretation of these functions 
is such that for a high coherence value, a value of the relative 
phase lead, or lag between the two signals, can be obtained. 

An overview shows, that at frequencies above and below 
the three frequencies studied where a strong possibility of a 


resonant bubble population existed, the phase and coherence 
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functions appeared much more confused. In all records, the 
Signals’ are in phase at low frequencies. For all the maximum 
phase shifts and zero crossings, the coherence values are very 
low. 

At the bubble frequencies of 56.2, 63.8, and 71.1 kHz, the 
coherence values are around 0.6 to 0.8 and the two signals are in 
phase within 30 degrees in the vicinity of the surface wave height 
frequency. This again underscores the affect of the sea surface 


wave height at these frequencies. 
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IX. SUMMARY AND CONCLUSIONS 


The statistics of the instantaneous output of two hydrophones 
separated by one meter were studied between the frequencies of 15 
to 151 kHz. The experiment was carried out from the NUC Oceano- 
graphic Research Tower off-shore at Mission Bay, San Diego. This 
study is of interest because bubbles and upper ocean fluctuations 
have Significant effects on the speed of sound in addition to 
those caused by temperature and salinity inhomogeneities. 

The experiment was conducted at a depth of 11 feet during sea 
state 1 conditions and a 4 ft peak to trough swell wave in the 
evening hours 19-2200 on 5 June 1972. The values ranged from 
+7.8 m/sec to -8 m/sec relative to the empirically predicted 
values of the Wilson and Leroy equations. The differential speed 
had a minimum corrected value of -2.4 m/sec at 19 kHz which 
suggests a batch of bubbles greater than 200 microns have been 
identified. Possible explanations are that these are bubbles that 
come from the bottom, or that they represent biological entities. 
There were speed of sound peaks at 41, 68 and 86 kHz. These 
identify a large population of bubbles, particuliarly in the 
range 56.3 to 71.1 kHz (bubble radii of 50 to 60 microns). 

The statistical analysis shows that for frequencies below the 
peak bubble populations the distribution of the instantaneous 
sound speed can be approximated by a Gaussian probability density 
function. At frequencies in and above the bubble populations, the 
variations in the sound speed lose the smooth Gaussian envelope 


and become irregular, sometimes double peaked. 
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The correlation functions show the surface wave height and 
phase fluctuations are in step at frequencies where bubble popu- 
lations appear. Above these frequencies the phase and wave height 
are relatively uncorrelated. There is evidence of long period 
correlation with internal waves. 

The spectral densities of both the phase and wave height show 
strongest values at frequencies of 0.4 Hz and below. The strongest 
values in the phase modulation spectrum were for sound in the range 


41 to 71 kHz and were identifiable with the surface wave spectrum. 
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